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We consider the self-evolution of strongly non-equilibrium interacting Bose gas. Due to the mere 
fact of large (as compared to unity) occupation numbers in the initial state the problem is directly 
reduced to the question of temporal evolution of the statistical matrix diagonal in the coherent-state 
representation. 

Strictly speaking, gauge symmetry is not destroyed even when the long-range ordering is com- 
pleted, owing to the ineviatble averaging over the ensemble. Actual symmetry breaking can occur 
only as a result of introducing a small term of the Hamiltonian violating conservation of particles, 
' or quantum-mechanical measurement, also implying non-conservation of particles. 



o 
o 



PACS numbers; 03.75.Fi, 32.80. Pj, 67.40.-w 



Kinetics of the Bose-Einstein condensate formation from a strongly non-equilibrium (completely disordered) initial 
1^ ■ state, ultimately leading to the long-range ordering, is one of the most interesting and fundamental dynamic problem 
C/^ ' of the physics of multi-particle systems. Recent exciting experiments on Bose-Einstein condensation (BEC) in ultra 
cold trapped atomic gases open up a unique opportunity of experimental study of this problem. Indeed, magnetic 
confinement completely isolates the system from the walls leading to the most refined situation, when the kinetics is 
the result of pure self-evolution of non-equilibrium gas. On the other hand, (forced) evaporative cooling technique 
^ I allows to have a variety of initial conditions differing by the number of particles and the residual energy. The formation 

• of the correlation properties can be traced non-destructively by measuring the secondary particles in the nonelastic 
[ processes, optical or dynamic response of the system. 

• Suppose that upon a fast cooling down there arises an isolated system with a certain fixed number of particles, 
I A'', the energy per particle, e*, being considerably less than the BEC temperature, Tc, that corresponds to the given 
. particle density, n. When such a system eventually comes to equilibrium in the course of self-evolution, the number 

of condensate particles will be close to N. The sharp contrast between this final state and condensateless initial one 
predetermines rather rich and nontrivial evolution scenario, involving a number of qualitatively different stages [^-Q. 
The particular analysis of the kinetics naturally starts with the statement that the evolution of the weakly interacting 
I . Bose gas with large (as compared to unity) occupation numbers is essentially classical-field phenomenon and thus can 
■ O ' be described from the very beginning in terms of non-linear Scrodinger equation. Taking this for granted, in Refs. 
^ [^-Q a detailed scenario of the BEC process was developed (see also Ref. |^). However, the very fact of possibility of 
the pure classical-field description is sometimes questioned. The heart of the doubts is the prejudice that for many- 
particle Bose systems the validity of the classical-field description implies the presence of the condensate. Recently, 
. 5^ ' a special effort has been undertaken by Stoof 1^ in order to substantiate the viewpoint that prior to the onset 
, of the classical-field regime there should take place a sort of essentially quantum phase transition - "nucleation of 
$H ' condensate" . (Published just now Ref. also contains the affirmation that there takes place a phase transition at 
an intermediate moment of evolution.) 

In this Letter we present a direct proof of the statement that the only fact of large (as compared to unity) occupation 
numbers in the completely disordered initial state forms the enough condition for reducing the problem to the question 
of temporal evolution of the statistical matrix diagonal in the coherent-state representation. This immediately leads 
to the substantiation of the classical-field description: At large occupation numbers the coherent states of the initial 
ensemble evolve following "independent trajectories" . Typical members of the ensemble are basically similar to each 
other, so for all practical purposes it is enough to consider the evolution of just one generic member of the ensemble. 
It is essential that though the eigenvalue of the field operator is non-zero for any particular member of the ensemble 
(and obeys non- linear Schrodinger equation), this by no means implies the presence of the condensate. 

Discussing the evolution of the correlation properties, we point out that the size of "long-range-ordered" region 
increases gradually with the suppression of the fluctuations of the phase, until the size of the system is reached. Thus 
we make sure that the condensation kinetics is a cross-over process, which does not imply any phase transition at 
some intermediate stage of evolution. 

Though one deals with an ensemble of states which are not eigenfunctions with respect to the gauge transformation 
of the field operator (virtual symmetry breaking), the genuine gauge symmetry is not destroyed, owing to the averaging 
over the ensemble. Symmetry breaking can occur only as a result of introducing a small term of the Hamiltonian 
violating conservation of particles, or quantum-mechanical measurement, also implying non-conservation of particles. 
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Let us turn to the formal analysis. For simplicity, we consider a homogeneous case. The parameter na^ (a is the 
scattering length) is assumed to be small enough to guarantee that the above-mentioned condition 



e, < Te (1) 

holds true along with the condition {h = 1) 

Ana , , 

nUo^e^ , Uo^ . (2 

m 

This means that at the initial moment most of the particles are in the kinetic region, where non-equilibrium state can 
be described in terms of distribution of occupation numbers. Clearly, for the characteristic occupation number, Uk, , 

{k1/2m — e^) we have 

..^(f^)"%>l. (3) 

To get a preliminary insight into the structure of the initial state, it is extremely instructive to estimate the density- 
density correlation function, K{r) — {h{Q)h{r)) , where n(r) is the number-density operator at the point r. As can be 
easily shown with the help of Wick's theorem, K{r) varies from 2n? (at r ^ 0) to at the distances of the order 
of A ~ 1/fc*. Together with Eq. (^) implying A 3> n~^^^ this means that the density experiences strong macroscopic 
fluctuations. Thus we realize that from the very beginning we deal with a sort of turbulent state with strong local 
non-uniformity. 

Further proceeding requires an explicit form of the statistical matrix. At the beginning, the statistical matrix of 
the whole system, W, is equal, to a good approximation, to the direct product of density matrices, T4^k, for each 
one-particle momentum state. The latter are diagonal in the occupation number representation: 

^ = , l^k = 5I^k(nk) |nk)(nk I . (4) 

k nk 

Since the coherent states 

|ak)=exp(--|akp) ^ \ n^) (5) 

form a complete set 1^, we can rewrite Eq. in the coherent-state representation. Assuming that Qk is a smooth 
quasi-continuous function of (a natural property of completely disordered initial state), one can prove that, up to 
smaller corrections in the parameter 1/nk, Wk is diagonal in new representation, satisfying the relation 



Wk = Jd{pl) J^Q^ipl)\a^){a^ 



(6) 



where ak — puexpUipy^). The most transparent way of proving this important result is to return from to (jij). Let 
us substitute Eq. (j^) into the right-hand side of the Eq. and perform the integration over ip\^. This leads to the 
expression diagonal in the occupation number representation: 

/ d{pl) QM E ^ /^k"" e""^ I "k ) ( nk I . (7) 

Then change the order of the summation and integration in the expression and notice, that the smoothness of Qk 
and large typical rik's allow replacing 

pl^^e-"'. ^ nu\S{pl-n^) , (8) 

immediately leading to Eq. (Q). 

Now we introduce the wavefunction of the coherent state of the system (at i = 0, a = {ak}) 

$(o,«) - ni"k> • (9) 

k 
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This function is an eigenfunction of the field operator ■0 with the eigenvalue 

^ = E/'ke'^-e'''"" . (10) 

k 

In terms of coherent states the statistical matrix of the initial state takes the form 

W{0) = J VaQ{a) \ $(0, a) ) ( $(0, a) | , (11) 

where 

Pa = n '^^pD — ' ^?(") = n Q^^pD ■ (12) 

k k 

Thus we see that at the beginning the statistical matrix of the system is diagonal in the coherent-state representation. 
Note, that Q is independent of (^Jk's. This means that the phases of harmonics in the initial state are completely 
non-correlated. 

Consider the evolution of the wavefunction <i>(t), which at the initial moment coincides with some coherent state 
$(0, a) After small temporal interval At we have 

$(At) = e"*-^'^*$(0) = $(0) - iAti7$(0) . (13) 

Here H is the Hamiltonian of the interacting Bose gas (the external potential V{t) is introduced for generality) 

H ^ ! drV3^[-— + V^(r)]v3 + — / dri^^^^p . (14) 
J 2to 2 J 

Let us act by the field operator on <I>(Ai). Since by definition i/>$(0) = V(0)«'(0), from Eq. (|l|) we have 

?/)$(At) = V'(0)«'(0) -iA#(0)iJ$(0) -iAi[Vi,7J]$(0) . (15) 

With the explicit form of the commutator, [tp,!!] = — (l/2m)A?/' + Vip + Uqi^'^i^i/j, the only nontrivial term in Eq. 
( p^ is -0^-0 -0$(O) = 'i/'^(0)'i/'^$(0). As can be checked straightforwardly, taking into account the smooth and rather 
narrow distribution of nk around Uk in the ak-state with = fik S> 1, we have 

V;t$(0) = V*(0)$(0) +C'(l/nk) ■ (16) 

So finally we obtain 

A 



V'(O) - lAt ( - — VXO) + V4>{0) + Uo I V(0) P V^O) 
' 2m 



$(At)+0(l/nk)+0((Ai)^) . (17) 



Repeating iteratively the above procedure and taking the limit At — > reveals that, apart from the 1/nk corrections, 
<i>(t) is a coherent state, 

= mm , (18) 

with eigenvalue ipit) obeying non- linear Schrodinger equation (NSE) 

= -^^^ + ViJ + Uo\i^\^^ . (19) 
ot 2m 

Hence, each separate coherent state of the ensemble evolves following an "independent trajectory" (cf. Ref. [^), 
and for the statistical matrix we thus have 

W{t) ^ j VaQ{a) \^{t,a)){^{t,a)\ . (20) 

Diagonal statistical matrix ( pO| ) with Q{a) > describes the temporal evolution of a statistical ensemble of coherent 
states $(t, a). Dealing with essentially macroscopic system, almost for all practical purposes it is enough to consider 
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the evolution of some typical representative of the ensemble with a particular choice of a = {ipu, Pu}- The solution of 
Eq. (|l^) for the initial state ( p^ ) defined by this set of parameters gives us the classical- field picture of self-evolution. It 
is worth noting, that calculating any correlation function with W{t) (|o|), one first performs averaging over a coherent 
state and only then averaging over the ensemble. Eventually, most of the gross-features prove to be insensitive to the 
initial parameters just after the first averaging. 

The average of the operator ip over the coherent state, 

(^) = ($(i,a)|V^|$(i,a)) = ijit,a) , (21) 

is non-zero and is defined by the solution of NSE (|l9|). However, ensemble averaging with W{t) (^) leads to 

{(i^)) = = . (22) 

Eq. (|l9|) and virtual broken symmetry for any component of the ensemble do not imply the presence of the condensate. 
Eq. ( p^ ) reflects global gauge symmetry preserving during all the evolution. 
Consider one-particle density matrix 

K{r,t) = i / dr'((^t(r + r')^(r'))) • (23) 



(fl is the volume of the system.) Averaging over a coherent state is straightforward: 

K^{r,t) - i y" dr'($(t,a) I V;t(r + r')^(r') I *(t,«)) ^ dr>*(r r', a) ^(r', t, a) . (24) 

To obtain this correlator, one should find the solution of NSE in the entire temporal interval. The analysis presented 
in Refs. has revealed that the evolution involves a number of stages. The process starts with an explosion-like 

wave in momentum space propagating towards the low energy region. At this first stage the potential energy is much 
less than kinetic one, and, taking advantage of the random-phase approximation, one can go over from the NSE to 
the kinetic equation, and to carry out rather detailed analysis, including numerical simulations 1^,0. The next stage 
takes place when most of the particles, which eventually will form the condensate, find themselves in the coherent 
energetic region e < Sc ~ noUo- In this region the potential energy becomes greater than the kinetic one, and NSE 
thus cannot be reduced to the Boltzmann equation. Apart from a short transitory period, the typical momentum 
of the coherent interval, ko{t), is much less than kc = v^mri<)I7o. As is shown in Ref. the solution of Eq. (|is|), 
parameterized as 

^ ^e'^ , (25) 
is characterized by small density fluctuations, ( (Sn/no)'^ ) ^ 1. Let us substitute ( p5| ) into 

K^{r,t) ~ y" dr'e-'x('"+""'^*'")+*x(''''*'") . (26) 

If at the moment t the characteristic width of the momentum interval is ko{t), then at r <C ro{t) = l/ko{t) the 
phase difference in ( p6| ) practically vanishes, yielding Ka ~ riQ. This result is insensitive to the ensemble averaging 
(p3|). Hence, at r < ro(t) the structure of density matrix corresponds to that of the condensate state. This so-called 
quasicondensate state does not, however, posses genuine long-range order: at r ^ '"o(i) the correlator (|2^) approaches 
zero. This follows naturally from the absence of phase correlations at such distances, which implies automatically 
anomalous character of long- wave correlations and the presence of irregular vortex structure Q . The absence of phase 
correlations means that the correlator 

(Vx(r)Vx(O)) « fc'IXkpe^*^- • (27) 

k<ko(t) 

vanishes at r ^ ro{t). Since the decay of correlations takes place at the scale ~ rQ{t), and is associated with the 
typical phase difference ~ tt, one estimates (7 ~ 1) 

IXk|'« ^ , k<ko{t) . (28) 
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Expanding phase difference in the exponent of Eq. (^6|) into the Fourier series, after standard calculations, we have 

= noe-^ , S = ^ fc' I Xk |' (1 - coskr) . (29) 

k<ko{t) 

According to (|2^), S can be estimated as 

S* « 47r7ln(fco(00 , fco(0?->l- (30) 

As fluctuations attenuate, transferring the excess of energy to the subsystem of normal excitations, fco monotonically 
decreases. When fco(t) reaches the minimal value fcmin ~ tt/X [L is the system size), variation of S with r stops and 
the correlator Ka takes on an independent of r constant value no. [Naturally, this result is not affected by averaging 
with the statistical matrix (^0|).] Now the long-range ordering is completed: we have the genuine condensate. 

As is discussed in Ref. [Q , the time necessary for this final stage of evolution depends on what particular process 
turns out to be the "bottle-neck" for the relaxation kinetics in the case under consideration. If this process is the 
relaxation of the vortex structure, then the time is ry oc L^. If the crucial process is the relaxation of non-equilibrium 
long- wave fluctuations of the phase, then the corresponding time is r^p (x (hydrodynamic regime), ox (x L 
(Knudsen regime) (see Ref. [Q for the details). 

Thus, we have demonstrated that the appearance of the condensate in the course of self-evolution of non-equilibrium 
system is an essentially cross-over process. The result ( p^ ) is preserved during all the time of evolution, indicating 
that the global phase of the condensate wavefunction is fundamentally unknown. The genuine broken gauge symmetry 
(and hence genuine complex order parameter) can appear at the final step of evolution only if the Hamiltonian contains 
a small term violating conservation of particles, or as a result of quantum measurement (which also implies exchange 
of particles in this or that form). 

This work was supported by INTAS (Grant No. INTAS-93-2834-ext), by NWO (Project No. NWO-047-003.036), 
and, partially, by the Russian Foundation for Basic Research (Grant No. 95-02-06191a). 



[1] M.H. Anderson et ai, Science 269, 198 (1995); C.C. Bradley et ai, Phys. Rev. Lett. 75, 1687 (1995); K.B. Davis et ai, 

Phys. Rev. Lett. 75, 3969 (1995). 
[2] B.V. Svistunov, J. Moscow Phys. Soc. 1, 373 (1991). 

[3] Yu. Kagan, B.V. Svistunov, and G.V. Shlyapnikov, Zh. Eksp. Teor. Fiz. 101, 528 (1992) [Sov. Phys. - JETP 75, 387 
(1992)]. 

[4] Yu. Kagan and B.V. Svistunov, Zh. Eksp. Theor. Fiz. 105, 353 (1994) [Sov. Phys. - JETP 78, 187 (1994)]. 

[5] Yu. Kagan, "Kinetics of Bose-Einstein Condensate Formation in an Interacting Bose Gas", Bose-Einstein Condensation 

(eds. A. Griffin, D.W. Snoke, and S. Stringari), p. 202, Cambridge University Press, 1995. 
[6] H.T.C. Stoof, Phys. Rev. Lett. 66, 3148 (1991); Phys. Rev. A 45, 8398 (1992); see also: "Condensate Formation in a Bose 

Gas", Bose-Einstein Condensation (eds. A. Griffin, D.W. Snoke, and S. Stringari), p. 226, Cambridge University Press, 

1995. 

Mathematically, the treatment of these Refs. is corrupted by an inadequate varying of the effective action - quite arbitrary 
rejecting the term with temporal derivative, which is responsible for crucial difference between the non-Euclidean action 
and a Euclidean one: Any uniform external-potential-like term in the non-Euclidean gauge-invariant action can be absorbed 
into the phase of the field, playing thus no role in the intrinsic evolution of the system. 

[7] H.T.C. Stoof, Phys. Rev. Lett. 78, 768 (1997). 

[8] R.J. Glauber, Phys. Rev. 131, 2766 (1963). 

[9] P. Carruthers and K.S. Dy, Phys. Rev. 147, 214 (1966). 
[10] D. V. Semikoz and L L Tkachev, Phys. Rev. D 55, 489 (1997). 



5 



